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Critical effects at complete and critical wetting in three dimensions are studied using a coupled 
effective Hamiltonian H[s{y),l{y)]. The model is constructed via a novel variational principle 
which ensures that the choice of collective coordinate s(y) near the wall is optimal. We highlight 
the importance of a new wetting parameter 0.{T) which has a strong influence on critical properties 
and allows the status of long-standing Monte-Carlo simulation controversies to be re-examined. 
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A suitable microscopic starting point for modelling 
fluid adsorption in systems with short ranged forces at a 
planar wall (situated at z = 0) is the Landau-Ginzburg- 
Wilson (LGW) Hamiltonian y 
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based on a local magnetization order parameter m(y, z). 
Here hi and c are the surface field and enhancement 
respectively. The bulk potential (/i(m) yields coexistence 
between bulk phases a and /3 (with corresponding 
equilibrium magnetizations ttIq > and < 0) for 
sub-critical temperatures T < Tc and zero bulk field, 
/i = 0. The magnetization at the wall is denoted by 
"^i(y) = ^(y, z = 0). The model is expected to 
show wetting behaviour so that for sufficiently large 
surface field hi{> 0) and high (sub-critical) temperatures 
T > T\Y (where Tw is the wetting temperature) the 
wall-/3 interface is completely wet by the a phase when 
h = 0~ . A simple surface phase diagram may then show 
two types of continuous wetting transition in which the 
mean thickness, {£), of adsorbed a phase diverges leading 
to large scale fluctuations in the position of the a/3 
interface, characterised by the perpendicular and parallel 
correlation lengths, £_± and ^||. We distinguish between 
critical wetting transitions which occur for T — > Tyy 
with fixed surface field (or for fixed temperature and 
hi — > h^~) and h ^ 0^ , and complete wetting in which 
{£) diverges for T>Tw (or hi > h^) and h ^ 0' . 

Beyond a mean field (MF) approximation, analysis 
based on (|^) has proved too difficult. Instead most au- 
thors use effective interfacial models based on a collective 
coordinate ^(y) which only accounts for fluctuations in 
the position (height) of the a/3 interface. The energy cost 
of the configuration is given by the Hamiltonian 



dy 



Sa/3(^) 



[V£(y)]2 + VK(£(y)) (2) 



where we have included the position dependence of the 
stiffness coefficient Eq,^(£) emphasized by Fisher and Jin 
(FJ) and not allowed for in the original capillary wave 



approach The direct interaction with the wall is 

modelled by the binding potential W{£) which is sensitive 
to the details of interatomic forces Renormalization 
group (RG) predict that for dimension d = 3 and short 
ranged forces the critical properties of both complete and 
critical wetting depend on the dimensionless parameter 
Lu(T) — ^g^"^ with K the inverse (true) correlation 
length of the bulk a phase and = T.ap{£ = oo), the 
free interfacial stiffness. The wetting parameter takes a 
universal value loc — 0.77 near bulk Tc §. As is well 
known, the predictions of strong non-universal criticality 
for critical wetting are not consistent with extensive 
Monte-Carlo work |^ which only show MF behaviour, 
and imply wfT) ~ 0. On the other hand more recent 
simulations of complete wetting reveal an effective 
value of uj(T) much larger than expected. 

Even within MF theory however the interfacial model 
(^ is problematic. While it correctly predicts cor- 
relations in the magnetization near (and along) the 
aP interface (in agreement with calculations based on 
the LGW model) it does not describe those occurring 
through the wetting layer i. e. from the wall to the 
a(3 interface. To remedy this Boulter and Parry (BP) 
1^,^ have introduced a second collective coordinate which 
models the relatively small magnetization fluctuations 
near the wall and results in a coupled effective Hamilto- 
nian. For complete wetting, occurring for T ^ Tyy, the 
nature of the second collective coordinate is unambiguous 
and may be efficiently modelled using an 'interfacial-like' 
field ^i(y) which accounts for (small) local translations 
of the magnetization profile. Remarkably for d — 3 
the coupling of fluctuations appears to have a profound 
effect on the critical behaviour and increases the effective 
value of to{T) in agreement with the latest Ising model 
simulations 0. This may be traced to the existence of 
a marginal operator in three dimensions corresponding 
to translations in the origin of the binding potential [0 . 
However, the BP analysis is not quite fully quantitative 
(even for T :s> Tw) and becomes pathological for 
temperatures close to critical wetting (T ~ TV) where 
the magnetization proflle becomes locally translationally 
invariant near the wall. The status of the second 
collective coordinate is then in doubt and the influence, 
if any, of the coupling of order parameter fluctuations 
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near the wall and a(3 interface cannot be assessed 
satisfactorily. 

In this paper we present details of a systematic and 
fully quantitative theory of three dimensional complete 
and critical wetting transitions, which reaffirms the 
important role played by coupling effects. There are 
two aspects of our work that we wish to emphasize. 
The first concerns the method by which the coupled 
effective Hamiltonian is derived from the underlying 
LGW model and relates to the choice of collective 
coordinate at the wall. Specifically we introduce a novel 
variational principle which ensures that the collective 
coordinate description most accurately approximates the 
functional integral over the phase space of magnetization 
configurations appearing in the partition function. We 
believe that this aspect of our work is of more general 
importance than the application to wetting described 
here. The second aspect of our analysis is the quanti- 
tative results for observables characterising critical and 
complete wetting transitions. These can now depend on a 
second dimensionless wetting parameter ^^^"^ , involving 
the stiffness coefficient of the collective coordinate near 
the wall. Within our theory the value of this parameter 
at Tw is explicitly determined as 
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with xii'' = ^bT^jI^ and are the surface suscepti- 
bility and (second moment) transverse correlation length 
of the surface layer at the wall-a interface, respectively. 
We estimate that for Ising-like systems fi(T) is close 
to unity (above the roughening temperature) and also 
approaches a universal value fl c as T — > Tc- Using 
reliable scaling and RG results [|l2t we predict that ^Ic — 
|wc yielding flc — 0.9. The influence of ri(T) allows 
a quantitative explanation of the hitherto anomalous 
simulations of Binder, Landau and co-workers 

Our method of allowing order parameter fluctuations 
at the wall is a perturbative one in which we first 
consider magnetization configurations which do not cater 
for such effects. To begin therefore we follow FJ Q and 
define the collective coordinate £{y) as the surface of 
fixed magnetization m'^ = 0. Next note that in order 
to construct the model Hamiltonian it is sufficient to 
consider the properties of planar constrained profiles |^] 
and we first focus on the profile mFj{z) (written m^(z; I) 
in [^) which minimizes HLGwinT] subject to the planar 
crossing criterion constraint 
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The second step concerns the inclusion of order param- 
eter fluctuations near the wall and for this we consider 
profiles which are similar to mp,j{z) but allow for local 
translations and enhancements of the magnetization. 



Note that these subtleties need not arise for the field 
(,{y) for which translations of the profile always dominate. 
Specifically we consider an additionally constrained 
profile to(z), defined as the minimum of (|l]) subject to 
both the FJ condition (0) and 



m{zi + ssin(5) = mFj{zi) ~\- KrUaSCOsS 



(5) 



where zi is an arbitrary position close to the wall (chosen 
such that < kzi <SC 1) which does not enter the final 
results. 
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FIG. 1. Detail of the planar magnetization profiles near 
the wall in scaled units. The broken curve shows 771(2) 
which incorporates a local enhancement and translation of 
the FJ profile (corresponding to the solid line). The proper 
coordinate s and angle 5 are shown. 

The simple geometrical meaning of the proper col- 
lective coordinate s and the variable S is illustrated in 
Fig. 1. Note that for <5 = f the fluctuation in the 
local magnetization represented by s corresponds to a 
translation of the FJ profile. The coordinate s(y) is then 
the 'interfacial-like' variable ^i(y) used by BP. However 
their description is clearly problematic if mFj{z) is 
flat near the wall (as it is close to Tw) and hence 
translationally invariant. To overcome this problem we 
instead chose a value 5 — 5* {I) which maximizes the 
linear response ^|j,_q of the FJ profile to a change in 
proper collective coordinate s. This may be written 



dm 
97 



Km„ cos 5 



s=0 



dmpj 
dz 
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which implies that S* satisfies 



dm 



(6) 
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This variational choice of angle S can be seen to maximize 
the 'difference' between m(z) and mpj{z) for a given 
s, i. e. in Fig. 1, s is 'directed' along the normal to 
mpj at z = zi. The collective coordinate near the 
wall then provides the most efficient way of exploring 
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the phase space of relevant magnetization configurations 
appearing in the partition function. One may compare 
with the 'worst' choice of S for which s is 'directed' 
along the tangent vector to mpjiz) at zi — variation 
of s then leaves the FJ profile unchanged (to first order). 
The same criterion can also be derived algebraically by 
considering the structure of correlation functions [ pT| . 
We suspect that similar ideas may be pertinent whenever 
a collective coordinate description is used to model low 
energy fluctuations. Considering our own problem this 
approach allows us to develop a coupled theory of both 
complete and critical wetting transitions which is not 
possible in the BP description. These phenomena are 
considered separately below. 

(A) Complete wetting: The essential observation here 
is that in the approach to a complete wetting transition 
{h 0" for T >Tw or hi > hY) the angle 6* is largely 
independent of £ and may be regarded as a function of T 
and hi. We find 
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where we have used cma{T) = h^ ■ Deep in the complete 
wetting regime (T Tw or hi :$> hY) S* approaches 
^ so that the proper collective coordinate is 'interface- 
like' and the optimum model is basically the same as 
the earlier BP approach. We can therefore anticipate 
that the effective value of the wetting parameter is 
renormalised in accordance with this simpler theory [^,^ . 
On the other hand as we decrease the surface field (or 
temperature) and consider complete wetting transitions 
closer to the critical wetting phase boundary the angle is 
small ((5* - T - Tw or S* ^ hi - h^) and the proper 
collective coordinate is 'spin-like' (i. e. an enhancement of 
the profile in a vertical direction in Fig. 1). The absence 
of any interfacial or translation component implies that 
any increment to lo{T), associated with the coupling of 
interfacial fluctuations, vanishes. 

The effective Hamiltonian is constructed from the 
constrained profile (with S = 6*) using well established 
techniques [pPJTl[|. To calculate the critical exponents 
and amplitudes it is sufficient to ignore the position 
dependence of the stiffness matrix elements, in which case 
the coupled model is 



+W{£- ssinS*) 
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where W{£) is basically the same function as that 
appearing in the CW model The parameter r is 



related to the transverse correlation length of the wall- 
a interface by ^^^^ = |iJll|- "^^^ momentum cut- 
off of the £{y) field is the same as that in the CW 
model (i. e. A2 ^ ttk Unlike the earlier coupled 

model which is problematic near the cross-over to critical 
wetting the cut-off Ai for the lower field is well-behaved 
at all sub-critical temperatures. In particular near the 
critical wetting phase boundary Ai is essentially identical 
to the cut-off of the LGW model (Aq say) since the 
proper collective coordinate is 'spin-like'. The influence 
of fluctuations on the critical singularities of the complete 
wetting transition is manifest in the critical amplitude 
il 



lim 
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The RG analysis of the new Hamiltonian is very similar 
to that for the simpler theory [|j and as before we find 
9 = 1 + ^ but with Lj, the renormalised value of the 
wetting parameter, given explicitly by 



uiT,hi)=Lj{T) + 
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where the higher order term is negative jl^. Thus the 

optimized theory yields a precise perturbation expression 

for the surface field or temperature dependence of the 

renormalised wetting parameter. It follows that the 

increment to w(T) vanishes quadratically as T — » T^^, a 

result which can not be derived from the BP model. This 

is illustrated in Fig. 2 and compares favourably with the 

available Ising model data (see Fig. 1 of |1^ ). Note that 

for T^y close to Tc the second term in ( |ll| ) reduces to 

/i —h}^ 1 ^ _ 
^ ^ so that ui has a universal expansion about 
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FIG. 2. Schematic illustration of the effective value of 
the wetting parameter near a critical wetting transition 
(occurring at hi = hY)- The coupling of fluctuations 
leads to an increment to the CW result u which then 
vanishes quadratically as hi — * i. e. for complete 

wetting transitions occurring close to the critical wetting 
phase boundary. 

(B) Critical wetting: In many aspects the optimal 
coupled model yields results for critical wetting similar to 
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the CW and FJ theories |,|. That is, if one ignores the 
position dependence of the stiffness coefficients, the cor- 
relation length critical exponent v\\ is non-universal and 
dependent on w [Q. Including the position dependence 
of the stiffness coefficients there is a possibility that the 
transition is fluctuation induced first-order very similar 
to the stiffness instability mechanism of FJ. However 
there is behaviour associated with coupling effects which 
is not present in either of these models. This is most 
simply seen by calculating the Ginzburg criterion for the 
surface susceptibility xi = and comparing the result 
with that of CW theory . Along the critical isotherm 
(T = Tw or hi = hY and h ^ Q~) xi should show cross- 
over from mean field-like behaviour (xi ^ |^|~^) to the 

1 

true asymptotic criticality [xi ^ 1^1 ^"^^ with ~ 3.7 
for Tw ^ Tc) when the correlation length ^|| is close to 
a value satisfying 

^(T^H^) {i Ml + (e||A2)^] + T+W - l} 

= l + 17(T^)ln[l + (C„aAo)2] (12) 

Again this shows the presence of two wetting parameters 
in contrast to a CW approach (l4| which is equivalent to 

= 0. The upshot is that the true critical region is much 
smaller than previously expected. If we assume that both 
and £,wa are of order the (unit) lattice spacing of 
an Ising model (which is a reasonable estimate for the 
existing simulation studies |^ ) then the cross-over region 
corresponds to ^|| ^ 60 (lattice spacings) which is an 
order of magnitude bigger than the CW result. Following 
[Q we estimate that for the Ising model simulations the 
bulk field h (written -j in [||) needs to be ^ 6 x 10^^ 
before non-MF behaviour is observed. This is much 
smaller than the range of fields studied and may well be 
beyond the limits of present day Monte-Carlo techniques. 

In summary, we have developed a precise quantitative 
theory of three dimensional wetting transitions, incorpo- 
rating coupling effects. Deep in the complete wetting 
regime (T 3> TV) our theory recovers the earlier model 
of BP but also allows a new discussion of cross-over 
effects from complete to critical wetting. In particular the 
increment to the effective value of the wetting parameter 
is predicted to decrease as T ^ (or hi — > h^^^) 
in agreement with simulations |^]. Applications to the 
critical wetting transition reveal a significant change 
to the Ginzburg criterion implying that MF behaviour 
for the susceptibility xi should remain valid until the 
transverse correlation length is extremely large. This 
appears to resolve long-standing problems associated 
with Ising model simulations Q which only find MF 
behaviour for this response function. 
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